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Abstract 



In this work, we study the Lie-point symmetries of Kepler-Ermakov systems presented 
by C Athornc in J. Phys. A24 (1991), L1385-L1389. We determine the forms of 
arbitrary function H(x,y) in order to find the members of this class possessing the 
sl(2, R) symmetry and a Lagrangian. We show that these systems are usual Ermakov 
systems with the frequency function depending on the dynamical variables. 

CO 

\q ■ 1 Introduction 

O . 

The basic idea behind the symmetry methods is to reduce the order of differential equations 
under considerations as much as possible so that the integration can be done easily [11]. 
It is well known that a dynamical system is integrable if a sufficient number of sufficiently 
simple (e.g. polynomial) independent first integrals can be found. The knowledge of a first 
integral amounts to the reduction of the order of the integration procedure by one. For sys- 
tems of N second order differential equations this order has to be counted as 2N. Ermakov 
systems are time-dependent dynamical systems [9]. They contain one arbitrary function 
of time, the so-called frequency function, and two arbitrary homogeneous functions of 
dynamical variables. A trend in the latest developments on the subject is to focus atten- 
tion on some special features of subclasses of generalized Ermakov systems in which the 
frequency function may depend on time, the dynamical variables and their derivatives [5, 
6, 8]. These subclasses may be tailored to suit some particular applications or special 
purposes. A central feature of Ermakov systems is their property of always having a first 
integral [4]. This invariant plays a central role in the linearization of Ermakov systems [2, 
6, 10]. A class of dynamical systems was presented by Athorne [1] which can be regarded 
as perturbations of the classical Kepler problem or of an autonomous Ermakov system. 
Kepler-Ermakov systems generalize the usual Ermakov systems while preserving the prop- 
erty of being amenable to linearization [1, 3, 6]. In the next section, we study the Lie 
point symmetries of Kepler-Ermakov systems. We find the Kepler-Ermakov-Lagrangian 
systems possessing the sl(2, M) symmetry. 
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2 Symmetry generators 

Systems of second order differential equations often appear in classical mechanics. So we 
use the notations common in that field [11]. We consider the generalized coordinates q a 
as dependent variables, t as the independent variable, and denote dq a /dt by q a . Then, 

if = w a {q\q\t) , a,i = l,...,N, (1) 

is the general form for a system of second order differential equations. The corresponding 
linear partial differential equation of first order is 

which admits 2N functionally independent solutions <p a = f a (q a ,q a ,t) that are first 
integrals of the system (1). So every solution of (1) can locally be written as q a = q a (ip a ,t), 
where d(p a /dt = 0. 

When the operator A is written in the form of (2), we can write the infinitesimal 
generator of the Lie point symmetry admitted by (1) as 

X = ^V)| + ^V)^- (3) 
Its extension is 

and the symmetry condition is 

[X,A]=AA, (5) 

where A is a non-constant factor, in general. By writing (5) explicitly and matching the 
terms on both sides, we obtain the following equations 

= = A, (6) 
\ dt H dt dt v ' 



If (7) is written in full, it is 

dw a 

+ q a q b q%c + 2q a q% c - q c q b V a bc + q% t - 2q h rf th - r& = 0, (8) 



where the subscripts b, c = 1,...,N, and t denote partial differentiations with respect to 
generalized coordinates and time. 

Kepler-Ermakov systems are defined [1] as the system of equations 

— H + — 

r 3 x 3 

y + w 2 (t)y = -^H + ^g(y/x), (9) 



x + w 2 (t)x = ~-^H + —f(y/x), 
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where H is a function of x, y and r = (x 2 + y 2 ) . Here / and g are arbitrary functions 
of the indicated arguments and w is an arbitrary function of time. If H, f and g are zero, 
we have a time-dependent harmonic oscillator which occurs in many practical contexts. 
In the case that H is taken to be zero we have generalized Ermakov systems. 

In order to calculate the Lie point symmetries of system (9) we consider the generator 
of the group of point transformations 

d d d 

X = £(x,y,t)— + rfi(x,y,t)— + rj 2 (x,y,t) — . (10) 

By inserting (9) into the symmetry condition (8) we obtain a polynomial equation in x 
and y. The coefficients of all monomials of the form x m y n must be identically zero. This 
yields the following system of partial differential equations satisfied by £, 771 and 772, 

£,xx £,yy £,xy 0> Vlxx 2^ x 0; Vlxy ~ £,ty 0; 

V2xy ~ t,tx = 0, r) 2 yy - 2£ ty = 0, Vlyy = V2xx = 0, (11) 

where we use subscripts to denote partial derivatives. The solutions of equations (11) can 
be written as 

£ = n{t)x + 5{t)y + a(t), 771 = ai(x, t)y + (3i(x, t), 

V2 = a 2 (y,t)x + (3 2 (y,t), (12) 
where the functions k, 5, a, a±, (3\, a 2 and (3 2 must satisfy the conditions 
oi\xx = a 2yy = 0, /3i xx -2k = 0, (3 2yy -25 = 0, 

C-ix — 5 = 0, a 2y — k = 0. (13) 
The solutions of these equations are 

ot\(x,t) = 5x + <f>i(t), a 2 (y,t) = ky + (f) 2 (t), 

(3 1 {x,t) = kx 2 + fo{t)x + <p 5 {t), p 2 (y,t) = 5y 2 + Mt)y + Mt)- (14) 

Substituting (12) and (14) into the symmetry condition (8) and setting the coefficients of 
the terms linear in x and y to zero we obtain 

fo(t) = (&-c 1 )/2, cf> 4 (t) = - c 2 )/2, k = 6 = 0. (15) 

We note that, up to here, /, g, H and w are arbitrary functions of their arguments. 
Inserting (15) into the rest of (8), we find that 

Cl = C 2 = 0, (pi = (j) 2 = (j>5 = $6 = 0, (16) 

(xH x + yH y + 2H)^+ a +Awwa + A&w 2 = 0. (17) 

The last equation implies that the function H(x,y) is not completely arbitrary but must 
be a function on the integral surface of the partial differential equation 

xH x + yH y + 2H = C (x 2 + y 2 ) 3/2 (18) 
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where C is a constant. As a special case, for example, if C is zero, H must be in the form 

H = -^l, (19) 

y 2 

where h is an arbitrary function of its argument and the function a(t) must be the solution 
of 

a +4wwa + A&w 2 = (20) 

which is the third order form of Ermakov-Pinney equation. Then equations (9) have the 
symmetry 

^/n<91..<91..<9 

X = a(t)- + -a W x- + -„(t)y- (21) 

which is also a symmetry of 

x y 

xy-yx ^g(y/x) + —f(y/x) = (22) 

y x 



obtained [4, 8] by eliminating w(t) and H terms in (9). If w ^ 0, equation (20) can be 



integrated once to obtain, aa — \o 2 + 2a 2 w 2 = c\ which can be reduced to the Pinney 



tion (20) can be solved easily if w = 0. Hence the system (9) with w = has the symmetry 



equation, p + w p = by the change of variable a = p (t). On the other hand, equa- 



X = ( Cl t 2 + c 2 t + c 3 ) | + ( Cl tx + fx) §- x + (c lt y + | y ) | (23) 



which splits into three components 



^ 9 d d d ^ d x d y d „ d , A . 

G ^'m +tX T x +tV 8-y a ^ t 8t + 2lT x + \^ G ^8i (24) 

which represent conformal and self-similar transformations and time translation. These 
generators satisfy the commutation relations 

[G 3 ,G 2 } = G 3 , [G 3 ,G 1 ] = 2G 2 , [G 2 ,G 1 ]=G 1 . (25) 

By redefining the generators as G 3 = a±, G 2 = C2/2, G\ = 173, where 

[0-1,0-2] = 2oi, [0-1,0-3] = o- 2 , [o- 2 , 0-3] = 2o- 3 , 

we see that the symmetry algebra is sl(2,R). Actually, it is the characteristic algebra of 
generalized Ermakov systems [5]. We have to note that equation (22) can be integrated 
and the result is 

-(xy — xy) 2 + / [uf(u) — u~ 3 g(u)] du = const (26) 

which is known as the Ermakov-Lewis invariant and (23) is also a symmetry of the 
Ermakov-Lewis invariant [4]. 
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If C is a nonzero constant, the function H(x,y) must be the solution of equation (18), 
that is 



H = - 



h{x/y) , C 



2^3/2 



and the function a(t) must be the solution of 



(27) 



a + -Ca = (, 
5 



(28) 



where w = and £ is a constant. This equation can be integrated easily and the result is 



a(t) = cie m + c 2 e 



-fit 



/3 2 ' 



(29) 



where (3 = 2iy^ and ci, c 2 are constants of integration. Substituting (29) into (21) we 
obtain the symmetries 



Btd I 3 8t( 9 d \ 



d_ d_ 

dx ^ dy ) ' 



3 /3 2 9t 



(30) 



with commutation relations 

1 5 

[Gi, G 3 ] = — Gi, [Gi, G 2 ] = -f3 G 3 , 



[G 3 , G 2 ] — 



(31) 



The Lie algebra is again sl(2,R) which can be shown easily by choosing G\ = 2o"i/- v /5/3, 
G 2 = 2as/^/5f3, G3 = 02/2/3. We conclude that Kepler-Ermakov systems have the Lie 
algebra sl(2,R) if the function H is in the form of (27). 
By using (27) in (9) we have the equations 



x = —x 



y = -y 



c hjx/y) 

5 y2 (^ x 2 _|_ y 2 ^ 3 / 2 

C h{x/y) 

5 y 2 [x 2 + y 2 f /2 



+ 



+ 



f(y/x) 

X 3 

g{y/x) 
y 3 



(32) 



possessing the sl(2, M) symmetry. These equations can be considered as the equations of 
motions of a particle with unit mass and can be obtained from a Lagrangian function, 



L=\{x 2 + f)-^{x 2 + y 2 )-\ 



f{y/x) , g(y/x) 



+ 



- *(ar,J/), 



(33) 



if 



ff(y/x) + x z g'{y/x) = 0, 



(34) 
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where the prime denotes the derivative with respect to the argument and 
xh(x/y) d^> yh{x/y) 



8x ~ y 2 ^2 + ^2)3/2' 



dy 



y 2 (x 2 + y 2)3/2' 



(35) 



so that the force is derivable from the potential function. The function ^f(x, y) is integrable 
if 



h(x/y) 



Cpy 2 
x 2 + y 2 



(36) 



where Co is an arbitrary constant. As a result we have the equations of motion 



-x 



y = -y 



Co 



c _ 

T ~ \ x 2 + ^2)5/2 

C Cp 

5 ( x 2 + y 2)5/2 



+ 



+ 



f{y/x) 



g{y/x) 
y 3 



(37) 



for Kepler-Ermakov systems possessing the sl(2, R) symmetry that are obtained from the 
Lagrangian 



C 



Co 



-3/2 _ 1 

~ 2 



f{y/x) g{y/x) 



x- 



y 



(38) 



with the condition (34). 

In polar coordinates the Lagrangian and equations of motion can be written as 

T 1 / -2 , 2q2\ C 2 Co G(9) 

where G(6) = sec 2 6f(taii9) + esc 2 ^(tan^), and 

.. - 2 G{6) C Co 

r — rti ^ h — r r = 0, 



d{r 2 9) G'{9) 



dt 



+ 



0. 



The Ermakov-Lewis invariant (26) comes from the integration of (41) and is 

(•tan 8 



-i /<tan v 

I=-(er 2 f + J [uf{u)-u- 3 g{u)] 



du. 



(39) 

(40) 
(41) 

(42) 



Very recently [7] , it was shown that the Ermakov-Lewis invariant is the result of a Noether 
dynamical symmetry for a class of Lagrangian-Ermakov systems. A Noether symmetry 
is a Lie point transformation that leaves the action functional invariant up to an additive 
constant [11]. It is also a Lie symmetry of the corresponding Euler-Lagrange equations. 
If 



(43) 
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is the generator of a Noether symmetry, then 



dq h 



-L 



(44) 



is a first integral that satisfies X</? = 0. For the Lagrangian in (38) a Noether symmetry 
generator is G3 and the first integral determined by (44) is the Hamiltonian. 

Finally we calculate the generator of Ermakov invariant for the Lagrangian (38). In 
general, if the Lagrangian L and a first integral 99 of a dynamical system are known one 
can determine the corresponding symmetry of dynamical character [11]. If 



X = C(g fc , q\ t)| + V a (q\ <i k ,t)^- a + ^(q\ Q k ,t)-^, 



(45) 



is the generator of a dynamical symmetry, then there exist a first integral (p such that 
8 2 L 



dq a dq' 



dip 
"dqb 



(46) 



holds. Conversely, if ip is a first integral, then rf — q a ^ determines a dynamical symmetry 
which is called a Cartan symmetry [11]. We can now use (46) to determine the dynamical 
symmetry corresponding to the Ermakov invariant (26) and then have 



X = (xy - yx) [ y-^- - x^- ) + 



dx 



dy 



y{xy - yx) + -^giy/x) - -^f(y/x) 
y L x 3 



d_ 

dx 



x{xy - yx) + ^g(y/x) - \f(y/x) 
y x 



d_ 

dy' 



(47) 



which is not a point symmetry [11] because it is not possible to add it a multiple of A and 
thereby to get rid of all derivatives in the coefficients of d x and d y . 



3 Conclusion 

In this work we concentrated on the Kepler-Ermakov systems, with w (t) = 0, possessing 
the sl(2,]R) symmetry. We found that these systems are usual Ermakov systems with the 
frequency function depending on the dynamical variables. We obtained the Lagrangian 
for these systems and calculated the Cartan symmetry corresponding to the Ermakov 
invariant. 
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